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Abstract 

We study classical integrable systems based on the Alekseev-Meinrenken dynamical 
r-matrices corresponding to automorphisms of self-dual Lie algebras, Q. We prove that 
these r-matrices are uniquely characterized by a non-degeneracy property and apply a 
construction due to Li and Xu to associate spin Calogero type models with them. The 
equation of motion of any model of this type is found to be a projection of the natural 
geodesic equation on a Lie group G with Lie algebra Q, and its phase space is interpreted 
as a Hamiltonian reduction of an open submanifold of the cotangent bundle T*G, using the 
symmetry arising from the adjoint action of G twisted by the underlying automorphism. 
This shows the integrability of the resulting systems and gives an algorithm to solve 
them. As illustrative examples we present new models built on the involutive diagram 
automorphisms of the real split and compact simple Lie algebras, and also explain that 
many further examples fit in the dynamical r-matrix framework. 
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1 Introduction 



The study of integrable many body systems initiated by Calogero ^ , Sutherland |2j and Moser 
jS] is popular since these systems have interesting applications in several branches of theoretical 
physics and involve fascinating mathematics. See, for example, the reviews in jH El EI • The 
generalized Calogero models associated to root systems by Olshanetsky and Perelomov [H], the 
spin Calogero models due to Gibbons and Hermsen and the Ruijsenaars-Schneider models 
[TU] are also very important. Many approaches were applied to prove the classical integrability 
of these models. The most general among these are the construction of a Lax representation 
for the equation of motion and the realization of the system as a 'projection' of some other 
system that is integrable 'obviously', for example the free geodesic motion on a Lie group or on 
a symmetric space (see jHl El E] and the review The result of Babelon and Viallet jT3| 
linking Liouville integrability to the form of the Poisson brackets (PBs) of the Lax matrix is 
also relevant for us, since the determination of the r-matrices encoding these PBs in Calogero 
models led eventually to the questions investigated in this paper. 

The r-matrices of the An type spinless Calogero models found in El E| depend on the 
coordinates of the particles and are not unitary. They were re-derived by Avan, Babelon and 
Billey in jUj using Hamiltonian reduction of an auxiliary non-integrable spin extension of the 
Calogero model. The definition of the non-integrable model relies on a 'quasi-Lax operator', 
whose PBs (of the type fj2.10|) below) are encoded by certain unitary dynamical r-matrices that 
yield the Calogero r-matrices upon reduction. The unitary r-matrices in question were later 
recognized ^Hl to be identical to the fundamental solutions of the so-called classical dynamical 
Yang-Baxter equation (CDYBE) that appeared first in studies of the WZNW conformal field 
theory fOl ■ Subsequently, the geometric meaning of the CDYBE and the classification of 
its solutions was investigated by Etingof and Varchenko |2I], and dynamical r-matrices have 
since turned up in several contexts. 

The 'wonderful but mysterious' calculations of [EI were understood on an abstract level by 
Li and Xu, who then proposed a method to associate an integrable spin Calogero type model to 
any spectral parameter dependent |221 or independent |2SI dynamical r- matrix defined in pT] . 
In our view the essential point in [22l|2SI is that a 'quasi-Lax operator' of a non-integrable spin 
Calogero type system can be directly associated with any solution of the CDYBE, and upon 
Hamiltonian reduction this yields a model admitting a Lax representation in the usual sense. 
Li and Xu also constructed new spin Calogero models by applying their method to the rational, 
trigonometric and elliptic solutions of the CDYBE given in [21]. (The spin variables can be 
rendered frozen only in the An case.) A companion construction of spin Ruijsenaars models 
was presented in [21]. These results (see also [23122]) are significant and deserve further study. 

The main purpose of the present paper is to characterize the integrable spin Calogero 
type models that may be associated with the remarkable classical dynamical r-matrices found 
by Alekseev and Meinrenken in [2Z]. These trigonometric solutions of the CDYBE are in 
correspondence with scalar product preserving automorphisms of self-dual Lie algebras. In 
particular, they contain r-matrices labeled by the automorphisms of the Dynkin diagrams of 
the simple Lie algebras, and it seems interesting to enquire about the corresponding integrable 
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systems. In our investigation we shall proceed in such a way to underline that the association 
of the spin Calogero type models to the solutions of the CDYBE does not actually rely on 
the machinery of Lie algebroids that features in fl^ , but rather it can be done in a very 
direct manner. This does not mean that the algebroid-groupoid technique is not important, 
since it plays a crucial role in the factorization algorithm worked out in fl^ EH EE] to solve the 
resulting models in the general (including the elliptic) case. Still, as these techniques are not 
widely known, it is worth stressing that in some respects they are not essential. 

In our case the algebroid-groupoid technique can also be bypassed since it will be demon- 
strated that the models associated with the Alekseev-Meinrenken r-matrices ()3.2|) can be viewed 
as projections of the natural geodesic system on a corresponding Lie group. We first show this 
at the level of the equation of motion, which immediately yields a solution algorithm. Then we 
prove that the spin Calogero phase space obtained from the dynamical r-matrix construction is 
identical to a reduced phase space coming from the cotangent bundle of the corresponding Lie 
group, which is reduced using the adjoint action twisted by the underlying automorphism. This 
generalizes Reshetikhin's derivation [2H] of the 'principal trigonometric' spin Calogero model 
by Hamiltonian reduction, which corresponds to the identity automorphism. Certain rational 
spin Calogero models will be derived in a similar way, too, generalizing results in ESI EOl E] • 

As illustrative examples, we shall present new spin Calogero models built on the non-trivial 
involutive diagram automorphisms of the simple Lie algebras. We shall also point out that many 
other models can be obtained using scalar product preserving automorphisms of reductive Lie 
algebras, including for instance certain models found earlier by Blom and Langmann and by 
Polychronakos by means of different methods. The degenerate complete integrability of the 
spin Calogero models corresponding to the identity automorphism was explained in PB| from 
the reduction viewpoint (without making explicit the relationship to dynamical r-matrices) . It 
should be possible to generalize the arguments of to our examples, and one may also try 
to quantize our systems by quantum Hamiltonian reduction. 

In the next section, the construction of generalized spin Calogero models based on solutions 
of the CDYBE is reviewed. Subsequently, this construction is applied to those r-matrices that 
are non-degenerate in a technical sense, defined in Section 3, which ensures that their equation 
of motion is equivalent to the Lax equation ()2.17|) . This holds for the rational r-matrices fl3.1|) 
and for the Alekseev-Meinrenken r-matrices (13. 2|) . In Section 3 the Lax equation ()2.17|1 is shown 
to be a projection of the geodesic equation on the underlying Lie algebra or its Lie group in 
correspondence with the rational and trigonometric cases. Section 4 clarifies the symplectic 
aspect of this projection. In Section 5 the spin Calogero Hamiltonians associated with the non- 
trivial involutive diagram automorphisms of the real split and compact simple Lie algebras are 
presented. After a brief summary, further examples and other issues are discussed in Section 
6. In Appendix A the r-matrices fl3.2|l are proved to be uniquely characterized by their non- 
degeneracy property, which also leads to a simple description of all quasi-triangular solutions 
of the CDYBE in the standard compact case. Appendix B contains some auxiliary material. 

Our main results are Proposition 2 in Section 3 and Proposition 3 in Section 4 concerning 
the projection method, and the uniqueness statement of Proposition A. 2 about the Alekseev- 
Meinrenken r-matrices together with its corollary dealing with the standard compact case. 
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2 Spin Calogero models from dynamical r-matrices 



We below recapitulate the construction of spin Calogero type models developed by Li and Xu 
1221 121] generalizing the earlier work of Avan, Babelon and Billey fr7\ ITH] . The essential point 
is contained in Proposition 1 below, which is actually very simple. For our present purpose, 
we consider only spectral parameter independent dynamical r-matrices defined on a self-dual 
Abelian subalgebra, /C, of a self-dual (real or complex) Lie algebra Q. 

The self-duality of Q means that it is equipped with a non-degenerate, invariant, symmetric 
bilinear form ( , ), which is real or complex valued depending on Q being real or complex. A 
subalgebra K. C Q is called self-dual if the restriction of ( , ) to /C remains non-degenerate. (In 
the physically most interesting cases ^ is a real Lie algebra and ( , ) has definite signature on 
/C.) The 'scalar product' ( , ) gives rise to the orthogonal decomposition 

^ = /c + /c^, x = x,c + x^ yxeg. (2.1) 

A dynamical r-matrix associated with /C C ^ is by definition a (smooth or holomorphic) map 

R:IC^ End(^) (2.2) 

on some open domain /C C /C, which is required to be antisymmetric with respect to ( , ), 
equivariant and subject to the CDYBE j2I|. Using the adjoint representation of Q, for an 
Abelian JC the (infinitesimal) equivariance property reads 

ad« o R{q) = R{q) o ad« Vg G /C, k E K,. (2.3) 

We introduce dual bases {T*} and {Tj} in /C, and let g* = {q,T'^) denote the corresponding 
coordinates on /C. Below we use the notation 

np np 

v.R = te—, {x,iVR)Y) = r{x,—Y) yx,Yeg. (2.4) 

The CDYBE ('with coupling z/') can be formulated as the condition 

E,{R, X,Y) = VX, Y eg, (2.5) 

where z/ is a constant and 

E,{R, X,Y):= —[X, Y] + [RX, RY] - R{[X, RY] + [RX, Y]) 

+ iVR)Y) + (Vy,i?)X - iVx,R)Y (2.6) 

In the 'triangular' case z/ = and in the 'quasi-triangular' case u = 1. 

Let us identify /C and Q with their dual spaces by means of the scalar product, extend the 
bases {T*} and {Tj} of /C = 1, . . . , dim(/C)) to dual bases {T"} and {Ti,} of Q {a,b = 
1, . . . , dim(^)), and define /^^ := ([T^, Th],T^). Then consider the phase space 

M:=T*]Cxg* c^}Cx}Cxgc^{{q,p,^)} (2.7) 
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equipped with the direct product of the natural Poisson brackets (PBs) on T*/C and on Q*. In 
coordinates, 

{q\p,} = 5], {ea,u = /:.ec. (2.8) 

Let us now introduce the function L : Ai ^ Q hj 

L: + (2.9) 

One may call L a quasi-Lax operator, because it satisfies PBs of the 'St Petersburg type' up to 
an 'anomalous term'. This is the content of the following proposition. 

Proposition 1 |j22|. 123] , For any dynamical r -matrix R on an Abelian, self-dual subalgebra 
K, d Q , the Q -valued function L \2.9jl on M. verifies the Poisson bracket relation 

{Li, L2} = [i?i2, Li + L2] - V5^i?i2, (2.10) 

where ^k: '■ {(1,P,0 ^ ^ic th^ evaluation map and R12 = {Ta, RTb)T"- (g) T^, Li = L ® 1, 
L2 = 1 ® L. 

Incidentally, the proof shows a stronger result. Namely, by assuming the antisymmetry and 
the equivariance properties one can calculate that 

{Li, L2} - ([i?i2, L, + L2] - V^^Ru) = (e, E,{R, T^ T^))T, ® T,, (2.11) 

which shows that the CDYBE ((231) is equivalent to the PB (jTTTl|l . 

The message from ()2.1()j) is that if the derivative term was absent, then the ^-invariant 
functions of L would provide a Poisson commuting family (as usual in classical integrable 
systems). Thus one should impose the first class constraints 

= 0, (2.12) 

and then the ^-invariant functions of L yield a commuting family of gauge invariant Hamilto- 
nians. This is a universal way whereby one obtains candidates for classical integrable systems 
out of dynamical r-matrices. (For particular r-matrices some other constraints could kill the 
derivative term in ()2.10|) . too.) As with constant r-matrices, not all these systems are necessar- 
ily Liouville integrable, and one may need to search for specific symplectic leaves in the reduced 
phase space to obtain interesting examples. But all these systems are at least 'Lax integrable', 
since the PBs of L with the Hamiltonian 

Hiq,p,0 = h{Liq,p,0) (2.13) 
for a ^-invariant function h on Q take the following quasi-Lax form: 

{L,H} = [RV,L]-{V^^R)V, (2.14) 
where V is the ^-valued function on the phase space defined by 

V(q,p,0 = {'^h){L{q,p,0) (2.15) 
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with the natural ^-valued 'gradient' of h. The crucial point is that ()2.14|) yields a Lax equation 
for L upon setting = 0. The most important Hamiltonian is 

ML) = 1(L,L), (2.16) 

and after imposing ()2.12|) it generates the evolution equation 

L=[RL,L]. (2.17) 

In the examples considered later the operator R{q) is zero on /C and [R{q) + |) maps K.^ 
to /C"*" in an invertible manner. Assuming that R{q) vanishes on /C, ()2.16p gives rise to 

Hiq,p,0 = 1{P,P) - l{^±,R\q)^±) - + (2-18) 

This Hamiltonian is of the spin Calogero type since R{q) is actually either a rational or hyper- 
bolic function of q. The fourth term in ()2.18p is a Casimir function and the third one disappears 
by ()2.12j) . The corresponding evolution equation can be obtained directly as 

V 

Q = {q,H} = p~ -^ic, 
p = {p,/J} = -(i?(g)ex,(Vi?)(g)U, 
L = {U,H} = [^,^p + R{qf^^]. (2.19) 

One sees from ()2.9|1 that if [R{q) + |) is invertible on /C-*-, then the last two equations become 
equivalent to (2.17) upon setting C,k = 0. We impose this constraint and are only interested 
in the time evolution of the quantities that are invariant under the gauge transformations 
operating as 

{q,p,U,L) ^ {q,p,e^^^e-^,e^Le-^), (2.20) 

where k is an arbitrary /C-valued function. The transformation rule of L follows from the 
integrated version of ()2.H|1 . where in the notation we assumed that ^ is a matrix Lie algebra. 
The time evolution of the gauge invariant functions of L is not changed if one replaces (2.17) 
by 

L=[RL-K,L] (2.21) 

for any /C-valued function k on the (constrained) phase space. This will be used in exhibiting 
the geometric meaning of the integrable systems covered by the above formalism. 

Natural analogues of Proposition 1 hold also for spectral parameter dependent as well as for 
spectral parameter independent dynamical r-matrices defined |25 on an arbitrary subalgebra 
K, d Q. Formula ()2.9p remains essentially the same for all cases, the derivative V^^ in ()2.10j) is 
replaced by with xil^P^O = [QjP] + ^k. for an arbitrary self-dual /C [221 • Direct verification 
is equally easy in any case as in ()2.11|) . The (spectral parameter dependent variant of the) basic 
formula (j2.1()j] first appeared in fZl^j, without referring to the construction of L out of R. 
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The universal character of the relationship between ()2.9|) , ()2.10|) and the CDYBE was explained 
in 122 1211 as part of a general theory based on Lie algebroids. To be more precise, |22| deals 
with spectral parameter dependent and [23 with spectral parameter independent dynamical r- 
matrices. After realizing ()2.1()j) . the construction of integrable systems by imposing constraints 
arises immediately and can be found in all the references mentioned. 



3 Interpretation as projection of geodesic motion 

We here show that the time development of the generalized spin Calogero model defined by the 
Lax equation ()2.17|1 can be obtained as a projection of the free geodesic motion on a connected 
Lie group G with Lie algebra Q if the r- matrix is quasi-triangular {u = 1), or on the Lie 
algebra Q if the r- matrix is triangular {u = 0). In special cases this statement is already known 
131231211, but these papers do not explain the relationship to dynamical r- matrices. Our result 
given by Proposition 2 below applies generally, its simple proof sheds light on the geometric 
origin of the Lax operators ()2.9j) and also directly leads to an algorithm to solve the model. 

Our assumptions on R are that it is compatible with the decomposition ()2.1|) . vanishes on /C 
(which can be achieved acting by a gauge transformation on the solutions of the CDYBE) and 
is non- degenerate in the sense that the restriction of [R{q) + |) to /C"*" is an invertible operator. 
We saw that these assumptions ensure the equivalence between ()2.17|1 and ()2.19j) with ^/c = 0. 
In the triangular case one has the r-matrices 

fl.K^E„d(e). , ^ I (3.1) 

where the inverse exists on a dense domain, e.g., if Q is semisimple and /C is a Cartan subalgebra. 
The quasi-triangular r-matrices satisfying our assumptions all have the form R = R^ with 

tf: C ^ EndlS). , » :^ { ° ^^^^^^^^ ^ ^^^^^^^^ _ » ^1 (3.2) 

where 6 is an automorphism of Q that preserves also the scalar product, /C lies in the fixpoint 
set of 9, and the inverse that occurs is well-defined for a non-empty open subset /C C /C. These 
quasi-triangular dynamical r-matrices are due to Alekseev and Meinrenken |2ZI (see also [Hlj). 
Their uniqueness under the above mentioned conditions is a new result proved in Appendix A. 
For convenience, we introduce the notation 

fl".W:=fi^«) + 5 = {^l_,-,,-...|^,)-, (3.3) 

To study the quasi-triangular case, we take a connected Lie group G with Lie algebra Q and 
suppose that 6 lifts to an automorphism B of G (or work locally where such a lift exists). We 
restrict the considerations to an open submanifold G G G of 'regular' elements g that admit a 
'generalized polar decomposition': 

g = Q-\p)e'^p-^ with peG, qet. (3.4) 
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A geodesic in G is a curve g(t) subject to the equation 

I {gitr'm) = 0, (3.5) 

where dot also denotes 'time' derivative. We are only interested in geodesies that lie in G. 

Proposition 2. Consider a curve of the product form 

git) = e-i(p(t))e'^Wp(t)-i (3.6) 

with smooth functions q(t) G /C and pit) G G. Define the 'spin' variable ^±(t) by 

^±it) := i?+(g(t))~^M^(t) G /C^ with M = Mk + M^:= p-^p. (3.7) 

Then the geodesic equation IS. ,5]) implies the same time development for the gauge invariant 
functions of q, q and ^± as does the Lax equation \2.11^ together with p = q- 

We start the proof by noticing that 

g-'g = P{q - (l - 9-'e-^^'')M)p-' = p{q - (l - e-'e-^''^U^)M^)p-\ (3.8) 

By using the definitions in ()3.3p and ()3.7p . this is equivalently written as 

g-'g = p{q - i?t(?)e±)p"' = pLp'' (3.9) 

in terms of the Lax operator 

Liq,q,U) = q-K{q)C±, (3.10) 

which has the same form as ()2.9j) with q = p and C,k = 0. Continuing with ()3.9|1 . ()3.5j) takes 
the form 

= f^igit)-'git)) = PiL - [L, M])p-\ (3.11) 
Since M = Mjc + R^iq)'^^±, this in turn is equivalent to 

L=[L,Rliq)'^^ + MK:]. (3.12) 

The 'geodesic Lax equation' ()3.12j) is to be compared with the 'spin Calogero Lax equation' 
()2.17j) . which now reads 

L = [R'iq)L,L] = [Rliq)L,L] = [L.RliqfU - \ci]. (3.13) 

The difference of the evolutional derivatives of L defined by ()3.12|) and by ()3.13|) is an infinites- 
imal gauge transformation in the sense of ()2.21|) . and thus the proposition follows. 

It is worth remarking that originally the gauge transformations ()2.20|) arose from the first 
class constraints ()2.12|) . but now we can also view them as due to the obvious ambiguity of the 
parametrization of g in ()3.4|1 . It is not difficult to show that this parametrization is applicable 
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in a neighbourhood of e'^° G G if and only if the r- matrix ()3.2|) is regular at go ^ In fact, for 
any such qq & JC the elements of G in a neighbourhood of e'^° can be uniquely written as^ 

g = e-^(e^)e^e-^ (3.14) 

where q varies in /C around go and rj varies in JC^ around zero. By combining this with Propo- 
sition 2 one can in principle solve the spin Calogero equation of motion (for gauge invariant 
quantities) as follows. First, take an initial value (go,Po = Qo, (C±)o)- Next, construct Lq out of 
these data according to ()3.10|) . Set po := e G G, and thus obtain from ()3.5|1 . ()3.9|1 the geodesic 

g(t) = e'?°e*^o. (3.15) 

Then, decompose this geodesic according to ()3.14|) in the form 

g^t) = e-\p{t))e''^'^p~\t) with p(t) = e''(*), ri{t)e}C^, (3.16) 

which is possible at least for small values of t. Finally, determine ^±{t) from ()3.7|1 and find also 
L{t) from ()3.9|) . This procedure provides the solution of the spin Calogero equation of motion 
(constrained by = 0) up to a time dependent gauge transformation of the type in fj2.2U|) . 
but this is enough since we are interested only in gauge invariant quantities. 

The non-degenerate rational r-matrices ()3.1|) are related to geodesic motion on the Lie 
algebra Q. To see this let us consider a curve X{t) ^ Q having the product form 

X{t)=p{t)q{t)p{t)-\ p{t)eG, q{t)et. (3.17) 

This gives rise to X = pLp^^, where -R(g) is given by ()3.1|) and 

L = q-R{q)^± with = R{q)-^M±, M = + = (3-18) 

Then the equation of motion ^-^r^ = translates into the 'geodesic Lax equation' 

L = [L, M] = [R{q)L - M^, L]. (3.19) 

This is gauge equivalent to the 'spin Calogero Lax equation' fl2.17|) associated with the 'principal 
rational r-matrix' 1)3.11) by the construction in Section 2. 

The solution algorithm sketched above must be equivalent to the one in for the overlap 
of the cases considered (which actually include only the systems ()5.4p ). Our derivation by the 
projection method was inspired by the seminal work of Olshanetsky and Perelomov ^T] and it 
appears much simpler to us than the derivation found in Section 6 of |23j . 

For a simple Lie group G and 6' = id our Proposition 2 is consistent with Reshetikhin's 
identification fIE\ of the 'principal trigonometric' spin Calogero model (with Hamiltonian of 
the form in ()5.4|) or ()5.8p below) as the result of factoring the natural Hamiltonian system on 
T*G by the adjoint action of G. Here G is the set regular elements in the usual sense, and 
g (restricted to an open Weyl alcove) together with its momentum p and the spin variable 
in }C^/e^ provide coordinates on the factor space T*G/G. Next we generalize this result to 
arbitrary automorphisms. 

^One sees this by calculating the derivative of the map m : /C x /C^ 9 {q, r/) i— > 8~^(e'')e'^e~'' e G at (go, 0). 
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4 Reduction of T*G by twisted conjugations 



We have seen that the spin Calogero equation of motion ()3.13|) results from the geodesic equation 
on the Lie group G. The geodesic motion corresponds to the Hamiltonian system {T*G, Q, H) 
with the natural symplectic form VL and 'kinetic energy' H associated with ( , ). In this section 
we demonstrate that the spin Calogero phase space together with its Poisson structure and 
Hamiltonian can be identified as a Hamiltonian reduction of {T*G, Q, H). 

As might be guessed from ()3.fi|l and the result of |23 for 9 = id, the reduction is by the 
cotangent lift of the 'twisted adjoint action' of G associated with B. We trivialize T*G by right 
translations and identify Q* with Q with the aid of the scalar product. This gives the model 

T*Gc^G^G = {{g,J)\geG, JeQ} (4.1) 

with 

n = d{J,{dg)g-'), H=^{J,J). (4.2) 

We then consider the 0-twisted adjoint action defined by 

Ad® : {g, J) ^ {e-\k)gk-\ Q-\k)je-\k-^)) VA; G G, (4.3) 

where our notation is adapted to matrix Lie algebras as usual. The equivariant momentum 
map generating this Hamiltonian action of G, ip : T*G Q* Q, is readily found 

^■.(g,J)^9{J)-g-'jg. (4.4) 

Our definition in ()4.3|1 may look slightly cumbersome as far as the usage of the inverse of B is 
concerned, but this is how it fits with the form of the r-matrix chosen in ()3.2|) . 

We are interested in the Marsden-Weinstein reduced phase spaces^ (T*G)^=^/Gfj,, where fj, G 
Q*nip(T*G) is chosen arbitrarily and C G is its isotropy group. We recall that {T*G)^=^/G^ 
carries a non-degenerate Poisson structure induced by the projection of the symplectic form f2. 
These reduced phase spaces together encode the structure of T*G/G, since they correspond to 
the symplectic leaves (or unions of leaves) in T*G/G in a one-to-one manner, simply by means 
of projecting {T*G)^=^ C T*G to the G-orbits passing through it. For a textbook, see e.g. 

Let O denote the coadjoint orbit of G through — /i. It will be convenient to realize 
(T*G)^=^/G^ by using the standard 'shifting trick' as 

{T*G)^=,/G, ^ {T*G X O)^^, /G, (4.5) 

where O is equipped with its own symplectic form, a;^, and ip is the momentum map for the 

~e 

diagonal action. Ad;. , of G on T*G x O: 

Ad® : T*G x03{g, J, ^ ie-\k)gk-\ Q-\k)JQ-\k-^), k^k'^), \fk G G. (4.6) 

^The factor spaces that appear here are not everywhere smooth manifolds, but possess some singularities. 
This complication could be resolved, e.g., by restriction to a dense and open smooth component (the principal 
isotropy type for the group action) or by invoking the theory of stratified symplectic spaces and singular 
reduction Kffij . In this paper we disregard these technicalities, since they do not affect our main point. 
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The symplectic form on T*G x O is given by 

n'^ = n + uj'^, (4.7) 

and the associated momentum map reads 

^{g,J,^) = 0(J)-g-^Jg + ^. (4.8) 

The identification ()4.5|) means that the reduced symplectic structure can be induced from fi'-^. 

We wish to show that the symplectic leaves of the spin Calogero phase space defined by 
the construction of j23 12^1 outhned in Section 2 fill certain dense subsets in the Marsden- 
Weinstein reduced phase spaces ()4.5j] . The appropriate symplectic leaves result by performing 
the construction of Section 2 in such a way that the spin variable is restricted there to the 
orbit O from the beginning. To avoid notational complications, here we make the simplifying 
assumption that G is a compact, connected simple Lie group. Denote by the set of fixpoints 
of the automorphism B of G. In the compact case the generic orbits in G under the B-twisted 
adjoint action, which we call twisted conjugacy classes, can be parametrized by a certain alcove 
in a maximal torus K of G® (see [HHI and references therein). It is also known jSHl that 
(1 — ^~^e~^*^9) is invertible on /C"*- for generic elements <? G /C, where /C is the Lie algebra of K. 

Proposition 3. Consider an automorphism Q of a compact, connected simple Lie group G 
and a Cartan subalgebra IC of the Lie algebra of G^. Let G C G be the Ad-Q-invariant open 
submanifold that consists of the Q-twisted regular conjugacy classes represented by e'^ with q E K,, 
where the operator {l—9~^e~^'^'')\lC^ is invertible andq i— > e'^ is injective on the open submanifold 
/C C /C. Take an arbitrary fi G il){T*G) and denote by O the coadjoint orbit of G through —fi. 
Then the Marsden-Weinstein reduced phase space {T*G X (see is identical to 

the phase space of the spin Calogero system defined in Section 2 by imposing the first class 
constraint / fFT^) on T*IC xO C M 

To prove Proposition 3, the key point is to solve the momentum map constraint = ()4.8|) 
by first bringing g E G to its 'diagonal' representative e"^ with q E JC. After doing this, the 
constraint ip = becomes 

9{J)-e-''"^'{J)+^ = 0, (4.9) 
which is easily seen to be equivalent to 

e^ = and e{J) = J,c - Ri{q)U, (4.10) 
where R^{q) is defined in ()3.3p . The resulting manifold, S C (T*G x 0^^_^ given by 

S:={{e'',J,c-0-'Rl{q)^^,U)\qe}C, G /C, UeOnJC^}, (4.11) 

is the gauge slice of a global partial gauge fixing for the action of G on (T*G x O)^^^. By 
inspecting this partial gauge fixing one sees that 

{T*G X O) . Q /G ~ S/K ~ /C X /C X (O n IC^)/K, (4.12) 
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where K is the Lie subgroup of G corresponding to /C and we used the identification S ~ 
/C X /C X (C n /C"*-) defined by the parameters (g, Jic,^±) on S. In terms of these variables the 
restriction of ()4.7|) to S takes the following form: 

n'^\s = d{J,c,dq)+uj''\onK^- (4-13) 

By setting p := J^,-, this becomes identical to the restriction of the symplectic form of T*/C x (9 
underlying the PBs ()2.8|1 . which is obtained by imposing the first class constraint ()2.12j) that 
generate the factorization^ by K. This proves Proposition 3. 

Notice that upon setting p = Jjc, 0{J) in the gauge S ()4.1H] yields precisely the con- 
strained quasi-Lax operator ()2.9j) . It follows that the Hamiltonian also reduces as required, 
since Hig,J) = \{J,J) = \{e{J),e{J)) . 

Equation ()4.4j) implies that il){T*G) consists of the coadjoint orbits of G through the elements 
of /C^ C Q* , and it is also clear that (T*G)^=^ is dense in (T*G)^=^ for any fj, G ip(T*G). It 
would be interesting to further elaborate the relationship between the reduced systems arising 
from T*G and from the entire T*G. 

For more general reductive groups, the spin Calogero model constructed in Section 2 is 
essentially (up to a discrete symmetry) identical to the factor of T*G by G, where G consists 
of the elements g in ()3.4|1 and G acts by twisted conjugations. It is an open problem to explore 
the possible connections between the global structure of the reduction of T*G and the spin 
Calogero models for compact as well as for non-compact groups. In studying this problem one 
has to take into account the non-conjugate Cartan subgroups of that exist in general and 
all 6-twisted conjugacy classes, whose description is rather complicated even for 6 = id [10] . 

It is easy to show that the rational spin Calogero model corresponding to the r-matrix ()3.1|) 
is a reduction of the geodesic system on T*0, where G Q contains the adjoint orbits of G 
through /C. For a compact simple Lie algebra Q with Cartan subalgebra /C this statement was 
proved in jSj, without mention of r-matrices. See also [201 E] for similar results about rational 
spin Calogero models associated with Cartan involutions. 

Let Ih be an inner automorphism of G, operating as Ih '■ G 3 g ^ hgh^^ with a fixed h E G, 
and consider the automorphism O := o If^. The corresponding twisted adjoint actions of G 
on T*G are related by the formula 

Adf = Ch~i o Ad^ oCh VA; G G. (4.14) 

Here Ch '■ T*G —>■ T*G is the symplectomorphism Ch '■ {g,J) ^ {hg, hJh~^), which also 
preserves the kinetic energy. It follows from ()4.14j) that the symmetry reductions of the geodesic 
system on T*G by the 0-twisted and the G-twisted adjoint actions lead to reduced Hamiltonian 
systems that are isomorphic by a map induced from Ch- By Proposition 3, this implies the 
isomorphic nature of the integrable spin Calogero models that can be obtained by applying the 
method of jSZlES] to 'input data' related by an inner automorphism of ^. In fact, one expects on 
general grounds that only different cosets in Aut(^)/Int(^) correspond to significantly different 

^The singularities of S/K can be avoided by restriction to a dense open subset of O in (|4.12(l . like in |22l l2Hj . 
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integrable systems, and this could be confirmed directly as well without using the Hamiltonian 
reduction of T*G. 



5 Examples based on involutive diagram automorphisms 
of real split and compact simple Lie algebras 

As an illustration, we here present the spin Calogero systems associated with a real split or 
a compact simple Lie algebra equipped with an involutive automorphism induced from the 
Dynkin diagram. The systems corresponding to the non-trivial diagram automorphisms of the 
An, Dn and Eq algebras have not been described before, while those associated with the trivial 
diagram automorphism appear also, e.g., in [221 123 12H] • To keep the presentation short, in 
all cases we write down only the r-matrix and the Hamiltonian after imposing the constraint 
^ic = 0. See eqs. (5.4), ()5.8|) and (5.14), ()5.17|1 for the Hamiltonians corresponding to the trivial 
and non-trivial automorphisms, respectively. 

First, we introduce some notations relying on standard results about Lie algebras. Let A be 
a simple complex Lie algebra with Killing form ( , ), Cartan subalgebra 7i, set of simple roots 
n = {f^k} and positive roots $+. Choose root vectors X±^^^ normahzed so that (X^^., X.^^.) = 1. 
Any symmetry r of the Dynkin diagram of A extends to an automorphism of as well as to 
an automorphism of A by the requirement t{X±^^) = for any ipk G H. The diagram 

automorphism r of .4 commutes with the Chevalley automorphism a defined by a{X±^p^) = 
—X:^^!^. We fix root vectors in such a way that 



We denote by Hr the real span of the elements T^^ := [X^^, X_^^]. 

With the above conventions, the split (normal) real form of A is the Lie algebra Ar defined 
by the real span of the base elements T^^, and X±^p. To describe the example provided by the 
data Q := Ar, /C := Hr, = id, we expand the 'spin' variable ^ E Ar as 



X. 



(t(X^) and {X^, X_^) = 1 



(5.1) 




(5.2) 



Equation ()3.2|) gives the standard r-matrix 




(5.3) 



where q G Ti-r is such that this function is well-defined. Upon imposing the constraint ^7^, 
the spin Calogero Hamiltonian ()2.18|) takes the form 








(5.4) 
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The compact real form Ac of A is the real span of the base elements iT^^. together with 

i 1 

^(X^ + and := 

For Q := Ac, IC := 'iHr and 6* = id using the expansion 



Y^:=—{X^ + X_^) and := —{X^ - X_^) G $+. (5.5) 



i = i^w. + 2^ {V^Y^ + g^^) , (5.6) 
we have 

In this case the spin Calogero variables are (ig, ip,C,) E IC x K. x /C-*-, and the Hamiltonian reads 

2 



1 1 . 7] + 



sm" ^ 

where q E Hr is such that sin 7^ 0. 

Now let us consider an involutive automorphism r of ^ induced from a Dynkin diagram 
automorphism as outlined above. Notice that r preserves the split and the compact real 
forms as well as the real Cartan subalgebra Hr- Denote by A^, Af, Af, and Hf the 
corresponding eigensubspaces of r with eigenvalues ±1. It is known that A'^ is a complex 
simple Lie algebra with Cartan subalgebra 7i~^, and A~ is an irreducible module of A'^ whose 
non-zero weights have multiplicity one jSHl- Furthermore, (H^, A^) and {iHr, A'^) are the split 
and compact real forms of (H^,A'^) with the associated real irreducible modules A~ and A~, 
respectively. Convenient bases of all these spaces can be obtained from the above considered 
('re-normalized') Chevalley basis of A by the 'folding' procedure. We recall this in Appendix 
B, and next summarize the result only. 

Let A be the set of roots of (7-^+,^+) with associated root vectors X+ for a E A. Let 
r denote the set of non-zero weights for {1-C^,A^) with corresponding weight vectors X^ for 
A G r. We can choose these elements so that 

(X+ X+J = l = (X,-,Xr,) and a(X+) = -X^,, a(X,-) = -Xl,. (5.9) 

The folding procedure ensures that these bases are such that A'^ is the real span of together 
with the X^ for a G A, and A~ is spanned by 7i~ and the X^ for A G F. Let A+ and F+ 
denote positive roots and weights. The folding procedure also ensures that A'^ is the real span 
of i7i+ together with basis vectors 

F+:=-^(X+ + X+J and Z+ : = i= (X+ - X+J Va G A+, (5.10) 

and A~ is the real span of iH~ together with basis vectors 

Y,- := ^{X^ + X:,) and Z^ := ^{X^ - XZ,) VA G F+. (5.11) 
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Now we consider the system associated with the data Q := Ar, JC := 9 := r. We write 
^ G Ar as 

e = e4+e«. + E^«^" +EC^A with ^t^K- (5.12) 
According to ()3.2p 

= ^ E + ^ E C tanh ^X- , (5.13) 

where this function is smooth for q G Ti.^. After imposing the constraint = 0, we obtain 
from ()2.18|) the Hamiltonian 



4 ^ sinh^^ gX^Wr'SH./ 4 ^ cosh^ ^ 

aeA+ ^""^ 2 AGr+ ^^^^^ 2 



(5.14) 



In the corresponding compact case, with data Q := Ac, /C := iH^ , 6 = t, we parametrize 
^ E Ac as 

e = i^^. + i^.,. + E + + E (^A ^a" + ^A ) • (5.15) 

«eA+ AGr+ 
For iq E }C C iH^ , we then obtain 

R^m = ^ E ^ (^^^« - - ^ E ^ (^A ^A - Ca • (5.16) 

aGA+ AGr+ 

After imposing the first class constraints = 0, the spin Calogero variables are (ig,ip, ^) G 
X /C X /C^, and the Hamiltonian reads 

a6A+ ^^^^ 2 Aer+ "-"^ 2 

We finish by listing the positive roots A4. and weights r+ that arise for the involutive 
diagram automorphisms of the classical Lie algebras. 

If ^ = -Dn+ij then A'^ = Bn and the module A~ is isomorphic to the defining representation 
of Bn- Therefore 

A_|_ = {cfc ci, em\^ ^ k < I < n, 1 < m < n }, r+ = {e^ | 1 < m < n }. (5.18) 
One may realize as the space of real diagonal matrixes, g, of the form 

g = diag(gi,...,g„,0,0,-g„,...,-gi) and : g g^. (5.19) 

If ^ = A2n~i, then A^ = Cn and 

A+ = {cfc ± eu 2em \ l < k < I <n, l<m<n}, r+ = {ck ± ei\l < k < I < n}. (5.20) 
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Now Ti.^ can be realized as the space of real diagonal matrixes written as 

q = diag(gi, . . . , g„, -Qn, -Qi) and : g g^. (5.21) 

For A = A^n one has A^ = Bn with A+ as in ()5.18|1 and 

r+ = {cfc ± ei, Cm, 2em |1<A;</<?T., 1 < m < n}. (5.22) 
Now 7i+ consists of the real diagonal matrices 

q = diag(gi, . . . , g„, 0, . . . , -qi) and Cm ■ q ^ q-m- (5.23) 

In the A = Eq case A'^ = F^^ and the weights of its representation on A~ can be found in 
jnni- Of course, one can also construct an integrable spin Calogero model by using the third 
order diagram automorphism of D^. 

If in the above we allow all variables to be complex, then the models given by ()5.4|) and ()5.8|) 
(or respectively by ()5.14j) and ()5.17|l ) yield the holomorphic spin Calogero model belonging to 
the complex Lie algebra A equipped with the corresponding automorphism. In other words, 
we obtained different real forms of complex spin Calogero models using the real forms Ar and 
Ac of A. For general considerations about real forms of complexified Hamiltonian systems with 
the spinless Calogero models as examples, we refer to [33]. 



6 Discussion 

One of the original motivations behind the theory of constant and dynamical classical r-matrices 
is the potential application to build integrable systems using the r-matrices. The aim of the 
present paper has been to contribute to this theory by characterizing the systems that corre- 
spond by the construction of [221123 to the non-degenerate dynamical r-matrices fl3.2|) . (See 
also the uniqueness statement of Proposition A. 2 in Appendix A.) We demonstrated that these 
spin Calogero type models are projections of the natural geodesic system on the underlying 
Lie group G, where the automorphism enters through the twisted adjoint action of G. The 
Hamiltonian reduction picture that we developed permits to understand the exact solvability 
of the reduced systems, in principle from any point of view ranging from the proof of complete 
integrability to the explicit construction of the solutions. 

We presented new integrable models built on the diagram automorphisms of the real split 
and compact simple Lie algebras. To illustrate the considerable range of applicability of the 
dynamical r-matrix construction, we now mention further examples that could be developed 
utilizing r-matrices of the form ()3.2j) . For instance, one could apply the construction to the 
Cartan involutions of the simple real Lie algebras. Cartan involutions are singled out since 
their use guarantees the definite character of the kinetic energy of the corresponding models. 
New models are expected to arise from those Cartan involutions that are outer automorphisms 
of non-split real Lie algebras. In fact, this holds for the real forms su*{2n) and so{2p+ 1, 2g + 1) 
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(p 7^ q) of the classical complex Lie algebras and for some real forms of i^g- A rich set of other 
examples can be derived by considering direct sum Lie algebras composed of > 1 identical 
copies of a self-dual Lie algebra A. The scalar product of A induces a scalar product on 

g := A®A®---®A (6.1) 

in such a way that one obtains a scalar product preserving automorphism 6* of ^ by the formula 

e-.g^g, 9{ui,U2,...,UN)-=iTUN,TUi,...,TUN^l), (6.2) 

where r is a scalar product preserving automorphism of A. In the simplest case r = id and 6 
is the cyclic permutation automorphism. We have inspected this case by taking an arbitrary 
simple Lie algebra for A and taking /C as the diagonal embedding of a Cartan subalgebra of A 
into g. Then it turned out that the dynamical r-matrix construction reproduces (for A = An) 
certain generalized spin Calogero models found earlier in [221 IBB] by different methods. These 
examples and their generalizations for affine Lie algebras will be further studied elsewhere. 

A very interesting question that remains to be investigated is whether all spin Calogero 
type models that may be associated with arbitrary solutions of CDYBE fl2.5|l and its spectral 
parameter dependent version can be understood as Hamiltonian reductions of suitable 'obvi- 
ously integrable' systems. The quantization of the spin Calogero models is another important 
problem that appears to be largely open. 

Finally, let us discuss the generalization of the construction of Section 2 to dynamical 
r-matrices based on non-Abelian subalgebras, which were introduced in [21] and appeared 
afterwards in various contexts (see e.g. ESI references therein). In fact, the main 
examples of such r-matrices are the 'extended versions' of the Alekseev-Meinrenken r-matrices 
()3.2|) given on the full fixpoint set, say go C g, of the automorphism 9 by the formula [2ZI 

f /(adg) on go 

<.:6«-End(e). (6.3) 

with the analytic function f{z) = | coth | — K Now, even if is non-Abelian, one can still 
use ()2.9j] to define a quasi-Lax operator L{q,p,^) on T*^o x ^* — x x ^5 but then the 
anomalous term in ()2.1()|1 becomes V^-Rg^t with xiQ^P^O = [q^p] + ^Go- This statement can 
be found in j2Sl and is also readily verified similarly to Proposition 1. Therefore, in the non- 
Abelian case one has to impose the constraint x = to get rid of the anomalous term. If 
a non-Abelian ^0 contains an appropriate Abelian /C, then x = represents a larger set of 
constraints than ^k. = ^ does, and x = also generates a larger gauge group. In the end, 
further supposing for simplicity that g is compact, it is possible to show that the reduced phase 
space of the 'non-Abelian construction' defined by imposing the first class constraint x = 
on T*^o X g* is the same as the one obtained from the corresponding 'Abelian construction' 
that operates by imposing ,^ac = on T*/C x g*. (If g is reductive, then /C must be a Cartan 
subalgebra of ^0 for ()H.2j) to be well-defined. A slight complication for a non-compact 
real reductive g arises since ^0 admits non-conjugate Cartan subalgebras in general.) By the 
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argument just sketched, we do not expect the r-matrices based on non-Abehan subalgebras to 
lead to interesting new integrable models. This observation motivated us in the first place to 
restrict ourselves throughout this paper to dynamical r-matrices based on Abelian subalgebras. 
Our argument is closely related to the mapping between different r-matrices defined by means of 
Dirac reduction of the space of dynamical variables as studied in ^3] . The connection between 
the Abelian and non- Abelian variants of the dynamical r-matrix construction of spin Calogero 
type models is elaborated in more detail in 



A On the non-degenerate Abelian dynamical r-matrices 

In this appendix we show that the r-matrices ()3.2|) found in j2ZI can be uniquely characterized 
among the solutions of the CDYBE (with coupling u = 1) hj requiring the additional properties 
(fXTjl and that ensure the equivalence of the Lax equation ()2.17|1 with the Hamilton 

equation ()2.19j) upon imposing the constraint ^k: = and setting p = q. The constrained 
Hamilton equation always implies the Lax equation, but is not necessarily equivalent to it in 
general. For compact Lie algebras the non-degeneracy condition ()A.7|) holds automatically, 
which allows us to obtain a rather strong result about the solutions of the CDYBE in this case. 

Before presenting our results we introduce a convenient notion of Cayley transformation. 
Let R G End(^) be an antisymmetric linear operator on a self-dual Lie algebra Q such that the 
operators R± := R± |idg are invertible. (This hold for all antisymmetric operators if ^ is a 
real Lie algebra with 'scalar product' ( , ) of definite signature.) We call 

C := R+RZ^ (A.l) 

the Cayley transform of R. We then also have the inverse Cayley transformation formula 

R = ^iC + l)iC-ir\ (A.2) 

Note that C G End(^) is an 'orthogonal' operator , i.e., 

{CX, CY) = {X, Y) , VX, F G g. (A.3) 

In fact, (jA.lj) and ()A.2|1 represent a bijective correspondence between the 'antisymmetric' op- 
erators R for which R± are invertible and the 'orthogonal' operators C that do not admit the 
eigenvalue (Up to some signs, this becomes the usual Cayley transformation for Q = so(3).) 

Let /C C ^ be an arbitrary (not necessarily Abelian) self-dual subalgebra with a non-empty 
open, connected domain /C C /C. Assuming that 

i?: /C ^ End(^), q^R{q) (A.4) 

is antisymmetric and R±{q) are invertible for all q E IC, we introduce the Cayley transform 
C : /C — *• End(^) by pointwise apphcation of ()A.1|) . One can then verify 
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Lemma A.l. The map R is IC-equivariant, iff its Cayley transform C is fC-equivariant, i.e., 

(V[T,,]i?) (g) = [adT, R{q)] ^ (V[r,,]C) (g) = [ad^, C{q)] (VT G /C, Vg G /C). (A.5) 

R solves the CDYBE \2. 5]) with v = \, iff C satisfies the equation 

Ciq)[X,Y] - [C{q)X,C{q)Y] + {C{q) - l){C{q)X, (VC)(g)F) 

-(V(c(,)y-y)^C)(g)X + (V(c(,)x-xkC)(g)r = (A.6) 

for all q E K, and X,Y ^ Q. 

Here, we are interested in antisymmetric r-matrices, R: JC ^ End(^), defined on an Abelian 
self-dual subalgebra K, <Z Q that are compatible with the decomposition Q = JC + K,-^ and satisfy 
the following non-degeneracy condition: 

i?_(g)|;c:i is invertible Vg G ^. (A.7) 

The compatibility really means that R{q)lC C /C since then R{q)K,-^ C /C"*" follows from the 
antisymmetry. For simplicity, we assume that the mapping R : IC End(^) is smooth or 
holomorphic for real or complex Q, respectively. Without further loss of generality, we can 
impose the auxiliary condition 

R{q)\K = 0, yqeL (A.8) 

This can be achieved, because (as is well known and easy to check) the restriction of the operator 
R{q) to /C is required by the CDYBE to yield an arbitrary closed 2-form, F, on /C by 

F(g) = Y,{T^, R{q)T,)dq' A dq^ (A.9) 

ij 

with a basis Tj of /C and corresponding coordinates g* of g G IC. If ()A.7jl and ()A.8j) hold, then 
we can use the Cayley transform (jA.ll) to prove the following result. 

Proposition A. 2. The (smooth or holomorphic) JC-equivariant solutions of the CDYBE h2. ,5]) 
with u = 1, on an Abelian self-dual subalgebra IC G Q, that are non- degenerate in the sense 
jyl. 7| ) and satisfy also iA.^} are precisely the maps of the form 

R{q)X = ^{9e'^-' + l){{9e'^'> -l)\ic±y' X, VX G /C^, g G /C, (A.IO) 

where 9 is a scalar product preserving automorphism of Q for which 9X = X for any X E )C 
and the inverse in iA.10\) exists. 

Proof. By the Cayley transformation formula (jA.lll the assumption R{q)\ic = translates into 
the equivalent equation C{q)\ic = — id/c, and our aim is to determine C{q)\ic±- To do that, first 
notice that in our case the /C-equivariance property ()A.5|) reads 

[adT,C(g)] = 0, VgG/C,VTG/C. (A.ll) 
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Now, let us turn our attention to the equation ()A.6|) for C. When X,Y G /C, then ()A.6|) 
becomes an obvious identity and gives no information. Next, consider the case when X G /C 
and Y G A^"*". Then we easily get 

{yxC){q)Y = [X,C{q)Y]. (A.12) 

We get the same equation when X G /C"*" and F G /C, up to change of letters. The most 
interesting case occurs when X, F G /C"*". Then, using the above two relations one finds that 
IBI) can be written as 



C[X, Y] - [CX, CY] + 2[CX, CY]^ = 0, (A.13) 

which is in fact equivalent to the following system of equations 

[CX, CY]^ = [X, F]^, [CX, CY]^ = C [X, Y]^ . (A.14) 

Now, let us introduce the function 

9^:IC^ GL(/C^), q ^ 9^{q) := {C{q)e-^'') (A.15) 

Equation ()A.12j) with (jA.lljl implies that this function satisfies the differential equation 

(Vx^±)(g) = 0, VgG/C,VXG/C. (A.16) 

Hence, 0±{q) is necessarily constant, Oj_{q) = 9j_ (Wq G IC). We extend this constant operator 
onto the whole Lie algebra Q by the formula 

e-.g^g, ^:=( Jf'^ (a.it) 

on JC . 

Then, one readily checks that C{q) satisfies equations ()A.3j) . (jA.llll . ()A.12|) and ()A.14jl if and 
only if the extended operator 6 ()A.17|1 is a scalar product preserving automorphism of the self- 
dual Lie algebra g. We have just demonstrated that under our conditions the Cayley transform 
()A.1|) of R{q) has the form 

- { Z f^-^^) 

where 9 G Aut(^) is a scalar product preserving automorphism with the property 9\!c = id^. 
So, by applying the inverse Cayley transformation ()A.2j) we obtain the Proposition. Q.E.D. 

Now we wish to point out some consequences. First note that if the self-dual, Abelian 
subalgebra K. C g is maximal as subalgebra, then R{q)]C C /C follows from the equivariance 
property. Thus in this case the /C^-component of any non-degenerate solution of the CDYBE 
(with coupling u = 1) has the form ()A.10|) . while its /C-component is furnished by a closed 
2-form F according ()A.9|) . This gives a characterization of the non-degenerate dynamical r- 
matrices on the Cartan subalgebras of the (real or complex) simple Lie algebras. If /C is a 
Cartan subalgebra of a complex simple g, then all automorphisms that fix /C pointwise have 
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the form 9 = e^*^" for a constant w G /C, and Proposition A. 2 reproduces a result of Etingof and 
Varchenko (see Theorem 3.1 in |2I])- More interestingly, notice that Proposition A. 2 can be 
applied to the classification of the r-matrices on Abelian subalgebras of compact Lie algebras, 
since in this case R± are automatically invertible for any antisymmetric R. Hence we obtain the 
following 'dynamical analogue' of the well-known uniqueness result jlH] on constant r-matrices 
of compact Lie algebras. 

Corollary A. 3. The JC-equivariant solutions of the CDYBE \2. 5]) with u = 1 on a Cartan 
subalgebra IC of a compact simple Lie algebra Q are precisely the maps q ^ R{q) for which 
R{(1)\k. corresponds to a closed 2-form by II A . .9)) and 

^(g)k^ = ^coth(^adg+,|^^) (A.19) 

with an arbitrary constant v E K,. The smoothness of R requires that ip{q + v) ^ 27riZ for any 
root if with respect to (the complexification of)K,. 

To obtain the corollary we used that the automorphisms of a compact simple Q whose 
fixpoint sets contain the Cartan subalgebra /C are of the form e'^^ with some v E IC. It can 
also be checked that from all (non-degenerate or not) solutions of the CDYBE that exist in 
the corresponding complex case precisely those giving ()A.19jl survive the restriction to the 
compact real form. 

For clarity, let us finally remind that in the CDYBE ()2.5j) the constant u can be replaced by 
P'^u if one replaces R{q) by (3R{f3q) with a constant (3^0- If 7^ 0, then for a complex Q the 
value of u can thus be scaled to 1, while for a real Q the two different possibilities u = 1 and 
u = i arise in principle. We have chosen u = 1 since this is required for Proposition 1. (The 
same choice is needed for applications in the WZNW model |12]-) Up to gauge transformations, 
for a simple Q all non-constant solutions with u that are known to us belong to u = 1. 
However, for a compact Q the constant solutions have z/ = i, and to our knowledge it has not 
been investigated if = i is excluded for non-constant solutions of the CDYBE or not. 



B Folding by an involutive diagram automorphism 

For convenience, we here summarize the 'folding procedure' referred to in Section 5. 

We denote by r the non-trivial involutive automorphism of the Dynkin diagram of ^ G 
{An, Dn, Eq} as well as its extension to A and to the root system $. Let us decompose the set 
of positive roots in the form 

$+ = SU^Ur(^), (B.l) 

where the elements of S are fixed by r, and the 2-point orbits of r in $+ are represented by \I'. 
Since r is a Cartan preserving involution of A, any choice of root vectors enjoys the property 
r(X<^) = c<^Xt-((p) with some constants subject to c^Ct-(^) = 1. If the root vectors are chosen 
as described in Section 5, then these constants satisfy also the relations c^^ = and c^c_<^ = 1. 
The first relation follows since r commutes with the Chevalley automorphism a, and the second 
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one is implied by 1 = (X^,X_(^) = (r(X^), r(X_^)). As a consequence, = 1 and = c^(^) 
for any (y9 G $. We can decompose S as S = S"*" U with 

S± := G $ + I r(<^) =^,c^ = ±1 }. (B.2) 

For any G $ denote by the restriction of ip to Ti^ . It is known (and is easy to check) that 
7^ \lip G $ and different r-orbits in $ yield different functionals on Ti^ upon restriction. 

In fact, a system of positive roots of introduced in Section 5, is provided by 

(B.3) 

(B.4) 

VypGH+, (B.5) 



X±^ V^GS-. (B.6) 

It is easy to see that the generators provided by this 'folding procedure' have the properties 
claimed in Section 5. 

In the most complicated example A = A2n, we have = {e^ — | 1 < A; < / < {2n + 1) }. 
With the usual conventions r : ^— — e2n+2-fc (VA; = 1, . . . , 2n + 1). Then we obtain 

^ = { Cfc - Ci, efc - e2„+2-«, em - e„+i |l<A;</<ra, l<m<n}, (B.7) 

E = E' = {era- e2n+2-m \l<m<n}. (B.8) 
The restriction of the above elements of and to 7i+ gives 

Cfc — Q = Cfc — e/, Cfe — e2n+2-« = Cfc + e/, im — e-ft+i = e^, — e2n+2-m = 2em,- (B.9) 

In the usual realization = skn+i, t ■ ^^ (-l)^~°+^-E2n+2-fe,2n+2-a wit h the ele- 

mentary matrices Ea^b for a, 6 = 1, . . . , 2n + 1. Then one readily verifies ()B.9|) by using 
that ea maps Q = diag((5i, • • • , (52n+i) ^ H to Qa ioi a = l,...,2n + 1 and maps 
q = diag(gi, . . . , g^, 0, —qn, ■ ■ ■ , —qi) G 'H^ to for A; = 1, ... , ra. The claimed properties of the 
root and weight vectors can also be checked on this model, but, of course, all the calculations 
can be done in a model independent manner as well. 
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A+ = I G ^ U S+ }, 
and the positive weights of {H'^,A^) are furnished by 

r+ = {(^ I G ^ u S" }. 

The root vectors used in Section 5 can be taken to be 
X+^ = i=(X±^ + r(X±^)) V^GV]/, 

and the weight vectors can be chosen as 



1 



X±^ = — (X±^-r(X±^)) Vy^GvI/, X. 
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